Non-equilibrium dynamics: Studies of reflection of Bose-Einstein condensates. 
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The study of the non-equilibrium dynamics in Bose-Einstein condensed gases has been dominated 
by the zero-temperature, mean field Gross-Pitaevskii formalism. Motivated by recent experiments on 
the reflection of condensates from silicon surfaces, we revisit the so-called classical field description of 
condensate dynamics, which incorporates the effects of quantum noise and can also be generalized to 
include thermal effects. The noise is included in a stochastic manner through the initial conditions. 
We show that the inclusion of such noise is important in the quantitative description of the recent 
reflection experiments. 
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I. INTRODUCTION 

The usefulness of the Gross-Pitaevskii equation 
(GPE)[1] in describing the dynamics of Bose-Einstein 
condensates (BECs) has been unreasonable in its effec- 
tiveness. At first glance, this non-linear Schrodinger 
equation describing the evolution of a purely classical 
field, corresponding to the order parameter, or wave 
function of the macroscopically occupied state, should 
only have any validity in regimes close to equilibrium at 
temperatures approaching absolute zero. Even in this 
regime, the apparent neglect of quantum fluctuations 
could in principle be a severe limitation. In practise, the 
GPE has been very effective in describing a wide array 
of experiments, many involving dynamical processes far 
from equilibrium. This success is rooted in the fact that 
the quasiparticle excitations of the Bose condensed sys- 
tem are also the collective excitations of the system^- 
As long as modes of excitation are highly occupied (mode 
occupancy large compared to unity), their evolution can 
be described within the GPE framework. Thus the GPE 
starting from a non-equilibrium configuration, to at least 
a certain extent, contains within it a description of the 
excited states and hence a representation of a thermal 
component. To make this more concrete and to indicate 
how quantum fluctuations can be included, we briefly re- 
view the classical field method and applications to BECs. 



II. THE CLASSICAL FIELD METHOD. 

Classical field methods, although suggested early in 
the development of quantum mechanics Q, were most 
extensively developed in the context of quantum laser 
physics [HHI]. The name is somewhat misleading in the 
sense that the methods were developed to account for ex- 
plicitly quantum effects such as quadrature squeezing Q 
and originates from the expression of the density opera- 
tor in a basis of coherent states using a classical "quasi- 
probability" distribution function. If we assume for the 
moment that such a distribution exists, one can show 



that the action of any quantum mechanical operator on 
the density operator is equivalent to that of a classical 
operator on the quasi-probability distribution function. 
Thus, the quantum mechanical problem, which is de- 
scribed by the evolution of the density operator as pre- 
scribed by a master equation, is transformed into a (clas- 
sical) partial differential equation for the evolution of the 
distribution function. 

Two main approaches, using different quasi-probability 
functions, have been used in the treatment of 
BECs. Specifically, these have used the positive-P 
representation @, B and, as we shall in this paper, the 
Wigner function [J 0, [l(J, [Hj], which was originally for- 
mulated as early as 1932. The P-function is used to 
express the density operator in the form 



p{t) = J d 2 a P{a,a* 



t)\a)(a\ 



(1) 



where \a) is a coherent state and the function P(a, a* ,t) 
is interpreted as a quasi-probability distribution. The 
"quasi" in the nomenclature is due to the fact that P is 
not necessarily positive-definite. The subclass of repre- 
sentations where it is are referred to as the positive-P 
representations . 



A. The truncated Wigner method. 

In this work we do not use the P-representation, but 
instead the (multi-mode) Wigner function[12j defined as 



W(aj,a*,t) 
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with characteristic function 
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where is the usual creation operator for an individual where 
bosonic atom in the mode j. 

The equation of motion for the density operator 



ih 



dm 

dt 



H(t),p(t) 



(4) 



can then be written in terms of the Wigner function 
through the use of certain operator equivalences. For 
example, 
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(9) 



Parameters for the mass, m and s-wave scattering length 
a = 2.9 nm are taken to be those for 23 Na. 

The quantum fluctuations arc introduced by combining 
the real particle field ip( r ) with a field of virtual particles 
X(r) to create a total field ^(r, 0) where the virtual field 
is defined as 



a)p{t) 



\£A W( aj ,a$ t t), (5) 



etc. 



To then proceed one considers, not the Wigner function 
evolution directly, which would be computationally in- 
tractable, but the stochastic differential equations which 
describe a single trajectory of a specific realisation of the 
system through phase space. Taking sufficient of these 
trajectories then allows the reconstruction of the Wigner 
function. To convert to the appropriate stochastic dif- 
ferential equations, it is necessary to truncate the full 
evolution, neglecting terms involving cubic derivatives. 
Fortunately this truncation can be justified provided that 
the total particle density is large compared to the mode 
density, which is the case in our simulations. In this man- 
ner a stochastic differential equation, of a form identical 
to the GPE, is obtained for each mode in the low energy 
subspace (below some cut-off) of the problem. 

More heuristically, one can think of the truncated 
Wigner method as simulating vacuum fluctuations 
through the addition of appropriate classical random 
fluctuations to the coherent field of the BEC's initial equi- 
librium state. To illustrate this let us explicitly consider 
the mode expansion of the wavefunction ^(r, i) given by 



1 M 

X(r) = ~i=^2 *»' ex P( ik J ' r )- 



(10) 



The complex amplitudes Xj have a Guassian distribution 
with (x*Xj) — \ $ij an d (XiXj) — 0, meaning that on av- 
erage each mode is populated by half a virtual particle, 
so that the total number of virtual particles is approxi- 
mately M/2. The truncated Wigner method is then valid 
provided the competing conditions, that the number of 
modes M is large enough such that the simulations are 
cut-off independent and small enough such that the num- 
ber of virtual particles is much less than the number of 
real particles, are satisfied. 

This yields, as does the full derivation using the trun- 
cated Wigner representation, the stochastic differential 
equation 



dt 



h 2 k 2 
— -( 

2m 



1 



dr e - ** 



[U cxt + U \V\ 2 ] * 
(11) 

for each mode within the low energy subspace, with J7 ex t 
now the total external potential. The dynamical evolu- 
tion of the BEC can now be obtained by solving these 
equations using a Fourth-order Runge-Kutta in the In- 
teraction Picture (RK4IP) algorithm [13J. 



1 M 



(6) 



where V is the real space volume and aj (t) now the am- 
plitude of the mode with wavevector kj with the wave- 
function normalised to the total number of atoms N. The 
mode space is taken to be spherical with the magnitude of 
the cut-off wavevector sufficient to prevent Fourier alias- 
ing. 

At times t < we imagine a BEC in equilibrium held 
in a magnetic trap centred at (—Ax, 0,0) with potential 
profile 



C/trap = y [<4{x + Ax) 2 + u 2 y 2 + lo 2 z 



2„21 



(J) 



such that the equilibrium mean field state is given by 
solving the three-dimensional GPE 
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III. REFLECTION FROM A BARRIER. 

Recent experiments have been performed to demon- 
strate quantum reflection of BECs from silicon 
surfaces [3, [Hj], showing their usefulness as atom mir- 
rors and traps. Significant disruption of the condensate 
was sometimes seen in the reflected condensate however, 
especially for dense clouds at low app roach velocity. Pre- 
vious work based upon the GPE[16| was able to explain 
much of this behaviour based upon disruption due to the 
mean field non-linearity, although the reflection probabil- 
ity and presence of an s-wave scattering halo could not 
be explained in this framework. Here we simulate two of 
the experiments using the truncated Wigner scheme. 

At t = the harmonic trap is displaced by Ax along 
the x-axis so that it is now centred atx = y = z = 0. 
This now accelerates the BEC towards an abrupt poten- 
tial step of positive height U parallel to the xy-plane at 
x = 0. The "impact velocity" of the cloud with the wall 
is thus v x ~ uj x Ax. If the kinetic energy of the atoms 
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at the barrier (ss ^ SL ) is much less than U then all the 

2 

atoms are reflected. If > U then there is finite trans- 
mission. 

The purpose of this paper is to simulate these recent 
experiments on quantum reflection from silicon surfaces. 
The surface is characterised by a Casimir-Polder poten- 
tial, which involves a rapid potential drop at the sur- 
face. This is impossible to model within the truncated 
Wigner framework, or indeed with the simple GPE in a 
mode representation, since the potential drop accelerates 
atoms to very large wavevectors which would exceed any 
practicable cut-off limited by computer memory. The 
Casimir-Polder potential is very abrupt however, with 
all spatial variation occurring over a length scale of less 
than a healing length. It can therefore be approximated 
by an abrupt potential drop. This of course does not ob- 
viate the problem of the acceleration of particles to high 
wavenumbers, however, instead of a sharp potential drop, 
we can employ an abrupt potential barrier which has the 
same transmission probability as the potential drop. Fur- 
thermore, previous theoretical workjlq]. using a different 
representation of the GPE, has shown the BEC dynamics 
of reflection to be qualitatively independent of the actual 
form of the barrier. 

Transmitted atoms are absorbed by the wall, which 
is simulated using an imaginary potential in the region 
x > such that 

C/imag = ~CX X > 0, (12) 

where C is a positive constant. This introduces a damp- 
ing term into the equations of motion for the mode am- 
plitudes in the region x > removing atoms. This imag- 
inary potential should only remove real atoms however, 
and not the noise of the virtual particles. To prevent the 
quantum noise being absorbed by the wall, the wall must 
"emit" virtual particles at an equivalent rate, increasing 
each mode amplitude aj by 

AV Vm ™ 

at each time step At. Here AT^ is the volume around 
one spatial grid point and the complex amplitude W v 
has the statistical properties (W u ) = 0, (W^W^) = 0, 

and (W*W„) = - A ^W <^. 



IV. SIMULATIONS 

In this section we present simulations corresponding to 
two recent quantum reflection experiments. The earlier 
experimental configuration [lij we refer to as A and the 
later [Hi as B. 



A. BEC A 

For BEC A, the total particle number N = 3 x 10 5 
and the trapping frequencies are uj x = 2ir x 3.3 rad s _1 , 
uj y = 2tt x 2.5 rad s _1 , and uj z = 2tt x 6.5 rad s _1 . These 
parameters yield a peak density of uq = 2.0 x 10 12 cm~ 3 . 

Initially we perform simulations without the addition 
of quantum noise to investigate the role of interferen- 
tial disruption in the reflection of the condensate from 
the potential barrier. This therefore corresponds to a 
simple GPE simulation. The positioning of the barrier 
relative to the initial trap centre (Ax) is chosen so that 
the characteristic velocity upon reflection is 1.2 mm s _1 . 
Constant density surface plots are shown in Fig. {TJ) for 
various times through the oscillation cycle back to refo- 
cusing of the BEC at the initial release point. 

The density of BEC A is rather low and hence the 
non-linearity in the GPE is weak for these simulations, 
corresponding to the weak mean field effects of the inter- 
atomic interactions. This leads to only very weak disrup- 
tion of the condensate upon reflection, especially for high 
approach velocities when the condensate has only limited 
time to interact with the standing matter wave grating 
formed upon reflection (Fig.fl]), panel 2). For character- 
istic velocities below about 2 mm s" 1 some disruption 
is seen however and vortex lines are produced in the re- 
sulting condensate as indicated in the final panel. This 
excitation of the internal structure of the cloud removes 
energy from the longitudinal motion of the condensate 
leading to a damping of the centre-of-mass motion [It]]. 
Since the excitation is due to the interference of the con- 
densate with the standing matter wave, we refer to it 
as interferential disruption. The observation of such dis- 
ruption is consistent with the ex per iments where "ex- 
cited and sometimes fragmented" [14f condensates were 
observed for small characteristic velocities of the cloud. 

Although even in the absence of quantum noise disrup- 
tion of the condensate occurs at low approach velocities, 
no scattering halo is produced and there is little deple- 
tion of the condensate. The depiction of the condensate 
comes from spontaneous processes. To simulate these we 
use the full truncated Wigner approach which is imple- 
mented by adding random (classical) fluctuations to the 
initial state as discussed in Sec. (II). 

At low approach velocities there is little difference be- 
tween the simulations with and without quantum noise. 
Little halo formation is seen and there is only weak 
depletion of the condensate. This is unsurprising as 
strong scattering halos were not observed in the earlier 
experiments By increasing the characteristic veloc- 
ity (increasing Ax) however, increased depletion and the 
formation of a halo can be induced. This is illustrated in 
Fig-© which shows the condensate depletion as a func- 
tion of time over one complete cycle for a characteristic 
approach velocity of 2.1 mm s _1 . The inset shows a sim- 
ulated absorption image, with optical access along the 
^-direction, at the end of the completed cycle. The s- 
wave scattering halo is clearly visible corresponding to 
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FIG. 1: Constant density surfaces for BEC A at various times. 
The characteristic velocity at the barrier is 1.2 mm s _1 . The 
circular arrow in the final panel indicates the sense of the 
excited vortex state. 



the distinct depletion of the condensate at this higher 
velocity. 




Time/ms 150 

FIG. 2: Depletion of BEC A for a characteristic approach 
velocity of 2.1 mm s _1 . The inset is a simulated absorption 
image after one complete oscillation showing the s-wave scat- 
tering halo. 



B. BEC B 

In the later experiments the BEC consisted of a higher 
total number of particles, N = 10 6 , with a more elon- 
gated configuration having trap frequencies of oj x — 
2-7T x 4.2 rad s -1 , ui v = 2n x 5.0 rad s -1 , and uj z — 2tt x 8.2 
rad s _1 , yielding an initial peak density of no = 5.2 x 10 12 
cm~ 3 . 

Since this BEC is denser and more elongated (again 
approaching the reflecting surface along its long princi- 
ple axis), it is much more severely distorted through the 
interferential disruption. This is seen in Fig. (|3]) which 
displays constant density surfaces of BEC B at vari- 
ous times through a full cycle in the absence of quan- 
tum noise. The characteristic velocity in this case is 2.1 
mm s _1 indicating that the disruption persists to signifi- 
cantly higher characteristic velocities. The central panel 
of Fig. ([3]| shows the formation of the standing wave. The 
final panel, at the end of the oscillation, shows a highly 
excited condensate, containing many vortices. 

In the later experiments using this higher density con- 
densate, a distinct scattering halo was observed together 
with a strong depletion of the condensate. This halo 
is absent in the simple GPE approach which was used 
to generate Fig.©. Since the formation of the s-wave 
halo is a spontaneous process, quantum noise must be 
included. 

Fig.(@J shows a simulated absorption image at the com- 
pletion of the reflection cycle (120 ms) for condensate B, 
with a characteristic approach velocity of 3 mm s _1 to 
match that of Fig. (4(e)) of Ref . [l5| . The barrier height 
was chosen so that the reflection probability in the sim- 
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FIG. 3: Constant density surfaces for BEC B at various times 
through the reflection cycle. The characteristic approach ve- 
locity is 2.1 mm s _1 . 



ulation was equivalent to that in the experiment. The 
scattering halo is very pronounced and the agreement 
between the simulated and experimental images excel- 
lent, clearly demonstrating the role of quantum noise in 
the reflection experiments. 




FIG. 4: Simulated absorption image for BEC B with a char- 
acteristic approach velocity of 3 mm s _1 . The barrier height is 
chosen so that the reflection probability in the simulation was 
equivalent to that in the experiment and the image should be 
compared with experimental image in Fig. (4(e)) of Ref . [l5| . 



V. CONCLUSIONS 

The utilisation of abrupt barriers, such as the Casimir- 
Polder potential of a silicon surface, as atom mirrors has 
already been experimentally demonstrated[3]. The ex- 
tension to more sophisticated mirror or guide geometries 
will require sophisticated simulation. In this paper we 
have developed such a fully three-dimensional technique 
within the truncated Wigner framework. We have noted 
the importance of the inclusion of quantum noise in sim- 
ulating spontaneous processes such as those leading to 
the formation of the s-wave scattering halo. 

It is interesting to note that there are two compet- 
ing effects leading to disruption or decoherence in the 
condensate after reflection. At low velocities, especially 
for higher density condensates, mean field effects lead to 
the interferential disruption of the condensate, whereas, 
at high characteristic approach velocities scattering into 
the (incoherent) s-wave halo has a significant depletion 
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effect upon the BEC. This suggests that applications of 
BEC reflection may be cleanest for moderate approach 
velocities which are above the threshold for interferential 
disruption, but below that for the formation of the scat- 
tering halo. An alternative might be to impose a static 
magnetic bias field in the vicinity of the reflecting surface 
to turn off the interactions and hence the interference via 
a Feshbach resonance. 
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